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We propose a two-path vortex interferometry experiment based on the Aharonov-Casher effect for 
detecting the non-Abelian nature of vortices in a chiral p-wave superconductor. The effect is based 
on observing vortex interference patterns upon enclosing a finite charge of externally controllable 
magnitude within the interference path. We predict that when the interfering vortices enclose an 
odd number of identical vortices in their path, the interference pattern disappears only for non- 
Abelian vortices. When pairing involves two distinct spin species, we derive the mutual statistics 
between half quantum and full quantum vortices and show that, remarkably, our predictions still 
hold for the situation of a full quantum vortex enclosing a half quantum vortex in its path. We 
discuss the experimentally relevant conditions under which these effects can be observed. 



Introduction. — It has been long understood theoreti- 
cally that topology in two spatial dimensions (2d) per- 
mits anyonic exchange statistics of identical particles in 
addition to the bosonic and fermionic possibilities [T]. 
Non-Abelian anyons, for which different exchange oper- 
ations may not commute, have attracted great atten- 
tion, particularly as building blocks of a topologically 
protected quantum computer Over the past sev- 

eral years, the quest for realizing non-Abelian anyons has 
been focused on the fractional quantum Hall (FQH) sys- 
tem (e.g. at filling factor 5/2) and the two-dimensional 
chiral p-wave superconductor (CpSC) [Sj- Theoretical 
proposals for detecting non-Abelian anyons in FQH sys- 
tems have exploited quantum interference , leading to 
active experimental study. In contrast, analogous explo- 
rations in the latter system have been scarce [51 . In this 
Letter we propose an interferometry experiment for di- 
rectly detecting non-Abelian anyons in a CpSC. 

In 5/2 FQH proposals, an Aharonov-Bohm interfer- 
ence pattern is expected when e/4 quasiparticles take two 
paths that enclose a region of flux and an even number 
of e/4 quasiparticles. Their non-Abelian statistics man- 
ifests itself as the effacing of this pattern when another 
e/4 quasiparticle is introduced in this region. Extending 
this idea to the setting of a CpSC poses a conceptual 
challenge. Non-Abelian anyons in a CpSC are neutral 
Majorana modes bound to vortex cores [6,, requiring an 
alternative means of interference. We show that the key 
to resolving this issue lies in the lesser known Aharonov- 
Casher (AC) effect [7 , which is an elegant 'dual' to the 
Aharonov-Bohm effect in that the roles of charge and flux 
are interchanged [8]. This is our flrst central result. 

The simplest incarnation of a CpSC, the 'one- 
component CpSC,' has been recently proposed in the su- 
perconducting state on the surface of a strong topological 
insulator induced by proximity to a conventional super- 
conductor , where locking of spin and orbital degrees of 
freedom effectively ensures that only a single component 
participates in the condensate. The situation is more 
complex in a 'two-component CpSC with an unconven- 




FIG. 1. (color online) Geometry of the proposed setup. A 
mesoscopic CpSC ring (yellow) is connected to leads (grey). A 
super-current density js is maintained in the sample. Due to 
the Magnus force, vortices (red) flow from bottom to top and 
choose one of two paths (L and R) to circumvent the central 
hole. A measurable voltage is induced across the leads by the 
vortex current J„ via the Josephson relation. Charge Q on 
the central island (blue) is controlled by gates (not shown) . A 
vortex may be stabilized at the hole, inducing currents around 
the hole edges shown by dashed arrows. 

tional pairing similar to that of superfluid "^He-A [TD] 
due to the existence of two distinct spin components. A 
prominent solid-state candidate in this class is the per- 
ovskite structured material Sr2Ru04 (SRO) JT]. The 
presence of both components doubles the number of Ma- 
jorana states in the core of a full quantum vortex (FQV) 
and renders their exchange statistics trivial. Fortunately, 
the order parameter allows for a topologically stable half 
quantum vortex (HQV) with a single Majorana bound 
state. A HQV carries a quantized flux of h/Ae and acts 
as a vortex in just one spin component. It has been ar- 
gued that certain bounded geometries may favor HQVs 
energetically [T^]. Recent tour de force experiments in 
Budakian's lab at UIUC have found strong evidence for 
the presence of h/Ae vortices bound to the inner radii 
of mesoscopic SRO rings [13]. Inspired by these experi- 
ments, we tailor our proposal to ring geometries and con- 
sider the effects of two spin components. It is likely that. 
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due to stability conditions, itinerant vortices are FQVs. 
We show that, quite remarkably, under some reasonable 
conditions much of our proposed interference physics is 
preserved even when FQVs go around a HQV stabilized 
around a mesoscopic hole. This is our second central 
result. 

Proposed setup and mechanism. — The geometry that 
we propose is shown in Fig. [T] An annulus formed by a 
mesoscopic CpSC sample is connected to leads. A super- 
current density, js, is maintained in the sample between 
the leads. Within the hole, an island supports variable 
charge, Q, whose magnitude is controlled by gating. An 
applied magnetic field can induce quantized flux through 
the hole. The interferometry is brought about by a con- 
spiracy of several factors as described below. 

The supercurrent induces a transverse force on vor- 
tices, which in its simplest Magnus form is proportional 
to js X $yz [T4|, where = h/e* is the flux quantum 
carried by the vortex and z is the unit vector normal to 
the superconducting plane. In a one-component CpSC, 
e* = 2e. In a two-component CpSC, e* — 2e (4e) for a 
FQV (HQV). This force causes vortices that enter from 
the bottom of the sample to move toward the top edge. 
These vortices can circumvent the central hole of the an- 
nulus by one of two paths, with amplitudes and tji 
and associated partial waves {iPl) s^nd \iPr) denoting the 
internal states of the vortices in the system. In princi- 
ple, vortices can be guided around the annulus by etching 
paths in the superconductor that would pin vortex tra- 
jectories. The associated vortex current is 

Jv (X Kip + liflp + 2\tL\\tR\Re [e^^(^L|Vi?)] , (1) 

with if — arg(t^ifl). The last term represents the inter- 
ference between the two paths. As vortices move across 
the sample, they generate a superconducting phase differ- 
ence normal to their motion. By the Josephson relation, 
a voltage drop K oc J„ develops along js. A measure- 
ment of Vs, or equivalently, of differential resistance for 
flxed current, therefore constitutes an interferometry of 
the two vortex paths. 

As a tuning parameter for the interference pattern, 
we propose the contribution to the phase ip stemming 
from the Aharonov-Casher effect, namely, the geometric 
phase picked up by the itinerant vortices upon encircling 
the charge Q on the central island. This non-local effect 
entails a magnetic moment fi moving in a static exter- 
nal electric field E. The adiabatic transport of such a 
moment around a closed loop causes it to accumulate a 
quantum phase [71[TS] (fiAC — /MxE-dl/(Sc^). A fluxhne 
carries a moment density eoc^^v along its length. When 
it encircles a total charge Q (of any distribution), it ac- 
quires a geometric phase (j)AC = ^vQ/fi = ^nQ/e* [l6] . 
It can be shown that the phase (j)Ac is truly non-local in 
that this form remains unaltered even in the case of an 
external charge placed at the center of a superconducting 
ring inspite of screening effects [71 117j . Consequently, as 



a function of charge Q, Jy in Eq. ([T]) and so the voltage 
Vs would exhibit oscillations with period e*. 

This picture changes drastically when a flux is ap- 
plied through the superconductor so as to nucleate a non- 
Abelian vortex on the hole in the center of the sample. 
Now the interference term also contains a braiding op- 
eration of the itinerant vortices around this vortex. As 
discussed in previous work [51 |T5], when the itinerant 
vortices are themselves non-Abelian, this braiding yields 
(V'ilV'-R) = 0; the interference term in Eq. ([T]) vanishes 
and Vs becomes independent of Q. Analogous to the 
FQH case, we propose the obliteration of the Aharonov- 
Casher oscillations as the clear signature of non-Abelian 
statistics of vortices. 

For a one-component CpSC, a vortex carries the small- 
est flux quantum and these considerations suffice. How- 
ever, a FQV in a two-component CpSC contains two Ma- 
jorana modes. One might therefore question whether 
this interferometry would work if the itinerant vortices 
are FQVs. Furthermore, the coupling between these two 
modes would split their energies and combine them into a 
regular fermion. We now consider this situation in more 
detail and show that there is an experimentally relevant 
window of parameters for which our proposed interfer- 
ometry does still work when FQVs go around a HQV. 

Majorana modes in a two-component CpSC. — The 
Hamiltonian for the two-component 2d CpSC is H = 

^ ^*-|- ^ , where h ~ (J^ ~ /i)/ is the single parti- 

cle hamiltonian {fi is the chemical potential and / is the 
identity matrix), and 

A - ^VAe^'^/^d^ - idy, {a- ■ d)a2}e'*/^ (2) 

is the pairing term. Here cr = (ci, cr2j era) are the Pauli 
matrices, d defines the direction of pairing in the space 
spanned by the spin-triplet states, $ is the phase of the 
order parameter, and va is a constant gap velocity. We 
represent the electronic operators as a Nambu spinor, 
— (V'j, V'|; V'ti V"^); where i/ij. creates an electron with 
spin a. Expanded in quasi-particle modes of energy E we 
have * = J2eXe^e, where Xe = iuEt,UEi,VEt,VEi) 
is the wavefunction of the Bogoliubov quasiparticle de- 
stroyed by r^;. Particle- hole symmetry dictates that 
Te = rl^; and thus UEa = V^eo-- energy, 
Fq = Fq = 7 would be a Majorana fermion. 

The Hamiltonian respects several more symmetries 
in the absence of external fields. It is invariant un- 
der a global SU(2) transformation = ^ q ^* ^ "^'i 

d = Rg^d', where S G SU(2) and Rs is the correspond- 
ing rotation S S0(3). The Hamiltonian is also invariant 
under a rotation by angle (f) around the z-axis, compen- 
sated by $ ^ + (j>, and a suitable rotation of d. Finally, 
there exists a Z2 symmetry under the mapping d — d 
and $ — ^ <i> -I- TT. 
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It is this last Z2 symmetry which ahows the presence 
of half-quantum vortices. When encircling a HQV, the 
phase of the order parameter changes by tt. The extra 
angular momentum needed to produce a single-valued 
Cooper pair wavefunction is provided by rotating d to 
— d. We consider a configuration where d rotates in an 
"easy" plane normal to a unit vector n, being the max- 
imal symmetry breaking that still allows a HQV. Ap- 



plying a magnetic field will explicitly break the SU(2) 
symmetry. We focus on the SU(2) symmetric case so as 
to freely choose this axis without affecting the result. 

As the vector field d in the presence of a FQV does 
not contain any winding, the two zero energy Majorana 
modes can be found by locally rotating d to the y-axis, for 
which the Hamiltonian decouples into spin components, 
see Eq. ([2|. The two Majorana modes are found to be 



E 



fSd 



V-'T 




(3) 



where Sd ~ g-«*o- n/2 -^^rjth n the unit normal to the plane 
of d and y, cos (5 = d • $ is the local phase of the order 
parameter in the absence of the vortex, and /(r) is an 
exponentially-decaying function of the distance from the 
vortex, determined by the vortex profile of A [TH] • 

For a HQV, when d rotates normal to n we may choose 
the y-axis at the intersection of the plane of rotation 
of d and xy-plane, so that fi — (sin/3,0,cos/3). Then 
the single Majorana mode of the HQV, 7', is also given 
by Eq. ([3]) in one of the spin components, with = 
e'"''^ y/^ which takes n z [TO]. 

Now consider a HQV on the hole in Fig. [T] and a FQV 
in the CpSC with a d-vector determined by the local d 
prescribed by the HQV. The situation relevant to our 
proposal is of transporting the FQV around the HQV 
adiabatically. The Majorana mode at the HQV experi- 
ences a relative change of phase $ by 27r as the FQV 
encircles it. Therefore, 7' — > —7'. At the FQV, d ro- 
tates by angle tt, i.e. (5 — > (5 + tt. When /? = (n||z), 
Sd exp (-if 0-3) Sd = -icr^Sd- Also, $ $ -|- tt so 
altogether the Majorana modes at the FQV transform as 

71^71; 72 -S- -72- 

In the general case (3 ^ 0, S —i{cr-h)S. This 
is where our specific choice for n pays off: we see that 
<T-ri = e-*^'^2/2^ggi/3<T2/2^ ^j^j^ g»/3-T2/2 i^gjj^g ^ real-valued 

rotation, so that by rotating to a Majorana basis 



71 ^ g^f a. 71 

72 / V 72 



we find that altogether 



71^71; 72 -72, 7 



(4) 



(5) 



Thus, the statistics in the new Majorana basis is the 
same as when 72 and 7' are braided, while 71 is unaf- 
fected. These are the same operations as those for braid- 
ing two HQVs binding Majoranas 72 and 7' [18]. Simi- 
larly, therefore, we find (V'LlV'i?) — 0- 

While the above corresponds to the SU(2) invariant 
situation, symmetry-breaking terms such as a Zeeman 



term or spin-orbit interactions would split the degen- 
eracy between the two Majorana modes on the FQV. 
For two such modes 71 and 72, having wave-functions 
Xi and X2 respectively, the associated coupling term as- 
sumes the form Hi — zA7i72. It should be noted that 
this term is invariant under the transformation intro- 
duced in Eq. Q, suggesting a larger symmetry in the 
low energy sector. Among the terms of the single parti- 
cle Hamiltonian that contribute to A, the most relevant 
is the Zeeman coupling |16| . hz — —gfiB^ ■ ''■/2, which 
explicitly breaks the SU(2) symmetry of the underlying 
Hamiltonian. In Nambu space, the coupling takes the 

^ff 1!t ) ■ -'^^ leading perturbation 

theory [16 we find (x2|-ffi|X2> = (xi|-H^i|xi) = 0, and 



form Hi 



X^-i{X2\Hi\xi) ^-9^lBB■d/2. 



(6) 



We see that when B _L d no coupling is generated. This 
is in fact an exact result: When d lies in the xy-plane, 
the effect of the Zeeman term is a relative shift of the 
chemical potential of the two Majorana modes. We there- 
fore get different functions fi,2{f) associated with the two 
spinors Xi,2j but no shift in energy. When the Zeeman 
term dominates the energetics, it will render d perpen- 
dicular to B, so that the two Majorana modes decouple. 
However, the combination of spin-orbit interaction and 
Zeeman term conspire to split the two Majorana modes. 

Experiment. — Piecing together the various compo- 
nents, we arrive at the following steps. The parentheses 
refer to a two-component CpSC. i) In the absence of an 
applied field through the hole in Fig. [l] drive a super- 
current js; ii) Vary the charge Q on the island enclosed 
by the hole. The differential resistance associated with 
the current ought to show AC oscillations as a function 
of Q due to vortex interference. The period of oscillation 
would be 2e (4e for itinerant HQVs); iii) Apply a field 
through the hole to nucleate a vortex (HQV). The inter- 
ference pattern should disappear; iv) Increase the field to 
nucleate a second vortex (FQV). The interference pattern 
(for itinerant FQVs) should reappear. If the interference 
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pattern is detected in ii) and does not disappear in step 
iii), the vortices either do not possess or cannot maintain 
their non-Abehan character over the time of the experi- 
ment. 

Discussion. — Our proposal requires several conditions 
to be met. For the AC effect, the vortices must be well 
separated from the charge on the island in Fig. [l] So, 
the dimension L of the hole must be greater than the 
magnetic penetration depth. Our proposal requires the 
coherent quantum motion of vortices. There is evidence 
that vortices can act as quantum objects [201 HI]- The 
biggest obstacle for such motion is dissipation in the vor- 
tex core. However, in the range of energies that we con- 
sider (< minigap) the core dynamics is effectively frozen. 
Another impediment lies in overcoming vortex pinning 
by impurities. Fortunately, superconductivity in SRO 
requires some of the cleanest crystals, evidenced by the 
observation of a fairly regular vortex lattice [22] . In addi- 
tion, methods have been developed using the concept of 
a flux flow transistor to create an "easy" channel for the 
flow of vortices, thus considerably reducing the critical 
current for the vortex flow. An alternative for achieving 
quantum behavior would be to replace Abrikosov vortices 
with vortices trapped in Josephson junction arrays as in 
the original experiment |20| . 

Some conditions must be met for any statistical in- 
teraction to be observed. First, the system must main- 
tain adiabaticity. Therefore, the rate for the motion of 
the vortices must be smaller than the smallest energy 
gap in the problem. This is the minigap that separates 
the Majorana zero mode in the vortex core and the next 
midgap state, Hva/L. So, the velocity of itinerant vor- 
tices must satisfy u <C wa- This is the most important 
condition for a one-component CpSC. Second, the Ma- 
jorana modes must remain coherent during the braiding 
operations. In a two-component CpSC, the splitting A 
of the two Majorana modes on the FQV results in de- 
coherence. So, in order to be sensitive to the braiding 
of Majorana modes in the presence of a HQV the braid- 
ing operation must run over a time shorter than nh/X, 
therefore XL/{Trh) < w <C wa- 

We now estimate the applicable range of parameters 
in SRO, where m ~ 10 electron's mass, hvA ~ 10-3 eVA 
and g « 2.5 [TT]. The vortex velocity is related through 
the Josephson relation to the voltage drop Vs across the 
superconductor: a vortex crossing the sample changes 
the superconducting phase by 27r. So, the voltage sen- 
sitivity Vs when iV„ vortices move across the sample is 
given by T4/7V„ « Trhv/{eL) < 10""^ volt for L « 1 mi- 
cron. The conditions for adiabaticity and coherence of 
FQV Majorana modes yield B • d <C 16 G, which needs 
only to be satisfied locally over a coherence length at 
the FQV. In fact, the orientation of the d-vector is still 
an unresolved issue; while a large enough Zeeman field 
along the z-axis ought to overcome spin-orbit interaction 
to render d _L B [131 US] , recent studies show that this 



orthogonality might even be the natural configuration in 
SRO ^24j . Additionally, the temperature must be smaller 
than the minigap not to excite other vortex core states 
which may lead to decoherence, i.e. T < hv/^/L ^ 1 mK. 
Other constraints on temperature might be imposed by 
the coherence condition discussed above. We believe that 
these conditions, while stringent, could still be realized 
in SRO. 

Another issue in SRO is the existence of multiple lay- 
ers. It is possible that this renders the exchange statistics 
of HQVs Abelian. (The non-Abelian nature of excita- 
tions in SRO has also recently been questioned in other 
ways [2S].) Our proposal can be used to see if this is 
indeed the case. Multiple layers may also result in deco- 
herence. The interlayer tunneling amplitude of Majorana 
fermions, t, depends on microscopic parameters (this 
is different from tunneling amplitude of Cooper pairs). 
Consequently, in the absence of experimental studies, it 
is not clear whether the limit of one independent Majo- 
rana mode per layer would hold or in fact the Majorana 
modes would completely hybridize to form a chiral low 
energy Majorana channel along the c-axis. If tunneling is 
present but small, the extra energy scale would introduce 
a decoherence time ^ h/t for our experiment [16j . 

Conclusion. — The proposal advocated in this Letter 
presents a smoking-gun signature of non-Abelian statis- 
tics of vortices in a superconductor. Our proposal hinges 
on vortex interferometry based on the AC effect; a sys- 
tematic experimental study of this effect in CpSC's would 
in itself be an important advancement. We have argued 
that under reasonable conditions in a two-componenent 
CpSC, not only HQVs but also FQVs encircling station- 
ary HQVs would exhibit non-Abelian statistics due to 
the Majorana modes in the vortex core. Our predicted 
signature of such statistics in CpSC's is the clear sup- 
pression of the AC interference pattern in our proposed 
set of experiments. 
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SUPPLEMENTAL MATERIAL 

Here we discuss some details of the Aharonov-Casher 
effect, our derivation for the Majorana zero modes in half 
quantum and full quantum vortices in a two-component 
chiral p-wave superconductor, the splitting of Majorana 
modes in the full quantum vortex, and the effects of mul- 
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tiple layers in SRO. 

THE AHARONOV-CASHER PHASE 

The Aharonov-Casher (AC) effect is closely related to 
the more well-known Aharonov-Bohm (AB) effect. They 
are essentially the same effect, viewed from two relatively 
different reference frames and are both special cases of 
the more general Berry's phase. 

In the AB effect the quantum mechanical wave func- 
tion of a charged particle acquires a geometric Berry's 
phase, 

^AB ^ jr A ■ dl, (7) 

where A is the vector potential and dl the element of the 
path of the particle. For a charge Q circling a region of 
flux $ — h/e* , we have (pAB = 2iTQ/e* . This is nonzero 
even if the flux is completely confined to an area that has 
no overlap with the particle's trajectory. 

In the AC effect, a moving fluxline acquires a Berry's 
phase as it circles a region of electric charge. This is the 
same AB phase as above but in the reference frame where 
the charge is stationary. Using the Lorentz transforma- 
tion of electromagnetism, the simplest form of the AC 
phase is [7] 

(I^Ac ^ <j) X E ■ dl, (8) 

when a particle carrying a magnetic moment fi closes a 
loop in an electric field E. 

Let us explicitly show that this is the same as the AB 
phase for any shape of fluxline and irrespective of the 
charge distribution, as long as the flux is confined. It 
was originally shown in Ref . [7| that for a point charge Q 
at ro and a moment fi at r 

E X /X = Qc^A^, (9) 

where A^ is the vector potential of the moment at rg, 
the position of the charge. A fluxline is a collection of 
infinitesimal magnetic moments dfi along its length (with 



so that it decouples into the two spin-components, and 
we get two independent zero energy states (see next sec- 
tion for details). Here, for simplicity we also assumed a 



a linear density dfi = eoc^^ds where ds is an element of 
length of the fluxline). Using Eq. Q in Eq. ^ we find 
quite generally 

<^AC = I hAA^/ds) ■ (-<fl) (10) 

=1/"* 

= 27rQ/e*. (12) 

We used the fact that § dA^ ■ {~dV) is nothing but the 
flux, d^, due to the moment d/i. as seen by the charge if 
it were to go around in an inverted loop, but in the same 
direction as the fluxline. The integration over the fluxline 
then gives the total flux carried by the fluxline, since the 
flux is confined completely inside the loop. This result 
does not depend on the distance between the fluxline and 
point charge. So, it simply generalizes to any geometric 
charge distribution as long as flux is confined and has no 
overlap with the charge distribution. 

MAJORANA MODES OF HQV AND FQV 
Full quantum vortex 

Let us take the d vector to be constant, pointing at 
some arbitrary direction. We perform a SU(2) transfor- 
mation 

This maps d — > d' = Rs<i, where Rs is the SO (3) rota- 
tion corresponding to S, since 

5^((T • d)a2S* ^S\cT- d)Scr2 ^ <t ■ Rgd a^. (14) 

We choose 5 to be a rotation that takes d to the y-axis, 

S' = e'3('^-fi) (15) 

where the unit vector n is orthogonal to the plane of d 
and y and 5 = cos~^(d • y). After this rotation, d — > 
d' — (0, 1, 0). The resulting Hamiltonian is now diagonal 
in spin indices. 



(16) 



I 

vanishing small vortex core so that is constant every- 
where but at the origin. Using the well known result for 
a HQV, we can write the two zero energy solutions for 



H' 



2m 



ivAe-''^/'^{d^ + idy)e-''^/^I 
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H' 



Xi 



^ ( fvt 

V2 \ rvt 



X2 



where = (1,0)"^, = (0, 1)-^ and 

1 hinr) \i < Imvl 



(17) 



(18) 



and K = y 2m|^ — imu^|. Rotating back to the origmal 
basis, the Majorana modes are given by Eq. (5) in the 



manuscript. 



Half quantum vortices 



Let us now consider a d vector that rotates by tt in 
a tilted plane of rotation, as shown in Fig. [5] We use 
the standard notation of the Euler angles, taking d to be 
along the X-axis, which sets it to be 



d(r) — (cosQ!COS7 — cos /3 sin a sin 7, cos 7 sin a + cos a cos /3 sin 7, sin /3 sin 7). (19) 
Choosing 7(r) = 9/2, in the polar coordinates r = {r,d), the d vector rotates by tt as we go around the origin. 



d(r, 9) ~ ( cos a cos - — cos /3 sin a sin - , cos - sin a + cos a cos /3 sin - , sin /3 sin - 



(20) 



The rotation 
takes d — )• d' = -Rsd in the xy-plane. So, 



„i§cr-N 



pie/2. 



(T • d')(T2 











(21) 



(22) 



Therefore, there is a vortex only in the spin-up compo- 
nent. Then, the zero energy solution is 



Xi 



}_ ( fvt 

V2 \ 



(23) 



with an associated Majorana mode 7'. If we take y-axis 
to be at the intersection of the plane of rotation of d 
and the xy-plane, then N = y, and the 7' is found from 
Eq. (5) in the paper in the spin-down component with 
S — e^^'^y/'^. 



THE SPLITTING OF MAJORANA MODES IN 
THE FQV 

Due to the presence of two spin components, two Majo- 
rana modes are present on a FQV. While these Majorana 
modes are protected against many types of local pertur- 
bations, they will split by the presence of perturbations 
that break spin rotation symmetry. These include the 
Zeeman splitting, hz = /icr • B, and spin-orbit interac- 
tions, hso = >^so • (cr X p). 



We first consider the case that there are no external 
fields and no spin orbit interactions. When the d vector 
is directed along the y-axis, the BdG Hamiltonian decou- 
ples into the two spin components. One can then solve 
independently the BdG equation in each spin component, 
and then rotate back to the d direction using a rotation 
S 



XI 



f*vt 



X2^U 



(24) 



where U is given by Eqs. ( 13 1 and (151 
terms Hz and Hso 



We now add the 



The key point is the transformation 
properties of these objects under rotations. 
Starting from Hz, 




FIG. 2. The Euler angles and the plane of rotation of d for 
the HQV. 
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Hz = 



/i(T • B 
-^^{a■ ■ B)^ 



-M(5'^cr • BS*)^ 



(25) 



Since S^cr ■ BS* = cr • -R5B, B transforms as a vector. We can therefore take d to lie along the y-axis, and B in an 
arbitrary direction. We then get 



(26) 



The spin-orbit term when Aso| |z was considered in Ref. 26 and found not to lift the degeneracy of the two Majorana 
modes. Here we show that, to the leading order in perturbation theory, this is true for general Xao and any direction 
of d. Using our solution in Eq. (24) we have 

Hso = 



Aso • (o- X p) 

-Xso ■ {cr X p)^ 



S-U^o • (cr X p)S 

-{S^Xso ■ (t X p)S)^ 



(27) 



Since S^Xgo ■ (f x p)^ = a ■ Rs{p x Xso), we can find 
the splitting in the basis where d is along the y-axis by 
rotating p and Xso by angle 6 around n. Therefore, 
in general the overlap {xi\Hso\X2) will contain terms 
involving components of the integral J {f*pf)dr. Not- 
ing that / and /* differ by a constant phase, we find 
/ {f*Vf)dY (X /(9r|/P)rfr — since / vanishes at infin- 
ity. We conclude 



{xi\Hso\X2)=0- 



(28) 



EFFECTS OF MULTIPLE LAYERS IN Sr2Ru04 

Non-Abelian statistics in the absence of inter-layer 
tunneling 

Sr2Ru04 is a layered material, and when a HQV is 
threaded through the system it goes through multiple 
layers. In the absence of tunneling between the layers, 
there will be one Majorana mode per layer bound to the 
vortex. Here we show that even in this limit the physics 
does not decouple into layers, but in fact results in an 
even-odd effect in the number of layers. 

We denote the layer index by £ = 1, . . . , A^. We start 
with a HQV going around a second HQV. For this case, 
the Majorana modes on each of the layers for the two vor- 
tices acquire a minus sign. This is generated by the fol- 
lowing transformation on the zero-energy Hilbert space, 



N 



lai- 



(29) 



To explore the effect of this transformation on the in- 
terference term, we should understand whether such 
transformations for different itinerant vortices will com- 
mute |4], i.e. whether \U' ,U"] = 0. It is easy to see 
that for I even they commute, while for i odd they do 
not commute, giving rise to the obliteration of the AC 
oscillations only in the latter case. 



In the case that a FQV goes around a HQV, we denote 
by 7^ and 72 the Majorana modes on the FQV, and by 72 
the one on the HQV. The transformation U' associated 
with the FQV encircling the HQV is given by 



N 

n 



72,£72,£ 



(30) 



which takes 7^ ^ — >■ I'n, but 72, £ — > — 72,^ and 72 ^ — >■ 
— 72£, being the same transformation (per layer) as in 



the main text. As the transformations ( p9| and (301 are 
the same, the results are the same: the obliteration of 
the AC oscillations will only occur for N odd. 

Finally, for a FQV encircling a FQV, the associated 
transformation is 



N 



^(Vi) = n 



ii..n2,niii2A 



(31) 



By relabeling the Majorana modes, the case of N lay- 
ers here is analogous to 2A^ layers for a HQV encircling 
a second HQV. Therefore, no effect on the interference 
pattern is expected in this case. 

To conclude, when either a HQV or a FQV encircles 
a HQV, we predict an even-odd effect in the number of 
layers: non-Abelian statistics will be present only in the 
case of an odd number of layers. In contrast, when a FQV 
encircles another FQV, the braiding statistics is abelian. 



Effects of inter-layer tunneling on the Majorana 
modes 



Inter-layer tunneling can be introduced between the 
Majorana states on a vortex (we consider a HQV for sim- 
plicity). We note that the tunneling of Majorana modes 
is a different physical phenomenon from that of Cooper 
pairs. For instance, when the sample has phase coher- 
ence along the z direction (i.e. in the absence of inter- 
layer supercurrents) and contains a single chiral domain. 
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the Majorana tunneling amplitude is exactly zero, even 
though Cooper pairs are free to move in all directions. 
This can be understood in two ways pT . 

First, since the z direction only appears in the kinetic 
energy, the momentum along the z-axis is a good quan- 
tum number. For a given pz the chemical potential is 
shifted as fi ^ fi — pi /2m. So, there is a degenerate flat 
band of zero energy Majorana states in the weak pairing 
regime, fi — pl/2m > [26]. Equivalently, from the solu- 
tion to the Majorana modes above, we can see that the 
overlap amplitude of two Majoranas on adjacent layers 
is proportional to sm{A(j)/2), where A(j) is the phase dif- 
ference between the two layers. Phase coherence means 
A(f) = and therefore all the tunneling amplitudes van- 
ish. 

This degeneracy is only protected by the additional 
symmetry in this idealized example and is lifted by phase 
fluctuations and other symmetry-breaking factors. But 
this example illustrates that the interlayer Majorana tun- 
neling is subtle and quite different from Cooper pair tun- 
neling. 

Let us now examine the spectrum of Majorana modes 
when there is a finite interlayer Majorana tunneling. The 
tunneling term acquires the general form 



H 



telele 



(32) 



We look for eigenstates of the form = J^iLi '"^7^' 
isfying [^^, F^] = ET^ . When N is odd there is always 
a zero energy solution to this equation, as is clear by 
the particle-hole symmetry of the Hamiltonian. Writ- 
ing the zero energy eigenvector of the Hamiltonian as 
(ui, . . . ,un)'^ , it is easy to see that U2n = (n e N) and 

n 

U2n+1 = U2n-l{t2n-l/t2n) = H {hn' -1 / hn') , (33) 

n' = l 

where ui is chosen by the requirement of normalization of 
the eigenvector. When all tunneling terms arc non-zero, 
the Majorana fermion is spread through the odd num- 
bered layers, avoiding the energy cost associated with 
nearest neighbors. When the number of layers is even, 
there are two such choices for states which then further 
hybridize and move away from zero energy. As we argued 
before, for an even number of layers we do not expect 
non-Abelian statistics to play a role in the AC interfer- 
ence. 

When the number of layers is odd, a second quantity 
which becomes of interest to us is the width of the result- 
ing band W . This width is set by the matrix elements 



ti, W (X maxiti. At zero temperature, the width of the 
band is of negligible consequence for a proper choice of 
the Majorana fermion operator. At finite temperatures, 
de-coherence results due to the finite band, with a typical 
time-scale which goes to infinity as — >■ or T — >■ 0. 
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